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ABSTRACT 

Recent observations of dipole anisotropies show that the velocity of the Local Group 
(uq) induced by the clustering of IRAS galaxies has an amplitude and direction similar 
to those of the velocity of Cosmic Microwave Background dipole anisotropy (ucmb), 
but the difference |ug — i’cmbI is still ~ 170 km/s, which is about 28% of |ucmb|- Here 
we consider the possibility that the origin of this difference comes from a hypothetical 
large-scale local void, with which we can account for the accelerating behavior of type 
la supernovae due to the spatial inhomogeneity of the Hubble constant without dark 
energies and derive the constraint to the model parameters of the local void. It is found 
as a result that the distance between the Local Group and the center of the void must be 
(10 - 20)h~^ Mpc, whose accurate value depends on the background model parameters. 

Subject headings: cosmic microwave background - cosmology: large scale structure of 
the universe - observations 


1. Introduction 

From the observed dipole component of temperature anisotropies in the cosmic microwave 
background radiation (GMB), our peculiar motion relative to the CMB frame has been determined 
because the dipole anisotropy is connected with the motion by a Doppler effect. The velocity of 
the Local Group relative to the CMB frame is found to be 

t;CMB = (-25.2, -545.4, 276.5) = 612.0 (-0.041, -0.89, 0.45) km s“\ (1) 

where the x-axis and 2 ;-axis are towards the Galactic center and the Galactic pole, respectively, and 
the total velocity |r7cMB| is 612 ± 22 in the direction of {I, b) = (268 ± 3, 27 ± 3) (Rowan-Robinson 
et al. 2000; Lawrence 1999). 

The peculiar motion is, on the other hand, caused gravitationally by the inhomogeneous matter 
distribution around us. Unitil now the latter velocity of the Local Group vq has been derived on 
the basis of linear gravitational instability using matter distributions that have been taken from 
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optical survey of galaxies, the IRAS galactic redshift survey and cluster surveys, for a long time 
(Yahil et al 1986; Meiksin & Davis 1986; Lahav 1987; Rowan-Robinson et al. 1990; Strauss et al. 
1992; Schmoldt et al. 1999; Branchini et al. 1999; Rowan-Robinson et al. 2000). The problem of 
the convergence of the velocity has also been discussed in many studies (Juszkiewicz, Vittorio and 
Wyse 1990; Scaramella, Vettolani & Zamorani 1994; Plionis & Kolokotronis 1998). 

The most recent determination of vq was given by Rowan-Robinson et al. (2000) using the 
Point Source Catalogue (PSCz) IRAS galactic redshift survey, which was analyzed out to a distance 
300/i“^ Mpc (Hq = 100/i km s“^ Mpc“^): 

Vq = u(rLG) = 697.7 uq km s“^, (2) 

where no = (—0.257, —0.811, 0.482) represents a directional unit vector. The above two velocities 
i^CMB and ug have similar amplitudes and directions, but, when inhomogeneities out to 300/i“^ 
Mpc are taken into account, the difference |ug — i^gmbI reaches ~ 170 km/s, which is about 28% 
of I'I^cmbI- The origin of this difference may come from a hypothetical object behind the Galactic 
center, inhomogeneous biasing or nonlinear correction to perturbation equations, as was indicated 
by Rowan-Robinson et al. (2000). 

In this paper we investigate the contribution of a hypothetical local void in which the inner 
Hubble constant (Hq) is larger than the outer Hubble constant (H^) and the radius of the boundary 
is ~ 200(/i^)“^ Mpc, where Hq = lOO/i^ km s“^ Mpc“^ (see Fig. 1). It is assumed that we are 
near the center of the inner region (e.g. Tomita (2000a,b)). 

//Fig.l // 

Here let us briefly describe the motivations for introducing the local void and the significance 
of its model parameters. Firstly, the observed values of the Hubble constant appear to be inhomo¬ 
geneous: the local value (Hq) in the neighborhood is about 72 km s“^ Mpc“^ (Sakai et al. 2000; 
Freedman et al. 2001), but the global values (H^) in remote regions are smaller than Hq, by factors 
which reach even 20%. The global values have been derived using gravitational leasing (Keeton 
& Kochanek 1997; Courbin et al. 1997; Fassnacht et al. 1999; Williams 2000; Tada & Futamase 
2001; Kochanek 2002a,b,c) and the Sunyaev-Zeldovich effect (Reese et al. 2000), although they 
depend on cosmological models and include some uncertainties. This inhomogeneity of the Hub¬ 
ble constant strongly suggests the existence of a local void (on scales with the redshift z < 0.1), 
because, if Hq > H^ locally, the central region must necessarily be underdense with respect to 
the total matter density. A local void on scales ~ 70 Mpc has once be studied by Zehavi et al. 
(1998) in terms of the Hubble bubble. At present, the existence of an underdense region on such 
scales is not clear from optical observations, athough it has so often been suggested on scales ~ 300 
Mpc (Marinoni et al 1999; Marzke et al. 1998; Folkes et al. 1999; Zucca et al. 1997), and so we 
should expect the situation that the distribution of luminous matter (such as galaxies and clus¬ 
ters) is nearly homogeneous, while that of dark matter is locally inhomogeneous. This situation. 



- 3 - 


which seems to be contradictory may be realized through the difference of the feedback system 
in hierarchical galaxy formation between the central (local) region and the outer (global) region 
(Tomita 2002), since in the low-density region, more galaxies are produced owing to less supernova 
explosions and photoionization, than in the high-density region (Kauffmann et al. 1999; Benson et 
al. 2001; Sommerville et al. 2001). 

Second, the accelerating behavior of high-redshift supernovae (SNIa) was found by two groups, 
the High-z SN Search Team (Schmidt et al. 1998; Riess et al. 1998, 1999) and the Supernova 
Cosmology Project (Perlmutter et al. 1999). Its origin has usually been explained by considering 
the role of nonzero cosmological constant A or dark energies in homogeneous cosmological models. 
Because the existence of nonzero A or dark energies is not evident theoretically, however, the search 
for other explanations may be significant. In the model with Hq > H^, on the other hand, the 
underdense region (the local void) plays the role of vacuum, and supernovae there show accelerating 
behaviors, even though we have no A and dark energies. In our recent papers (Tomita 2001a,b), 
it was shown that, for larger {Hq — the magnitude (m) and redshift (z) relation of high- 

redshift supernovae can be well reproduced in the models with smaller A, and, if {Hq — Hj^)/H^ = 
0.2, the relation can be quantitatively explained in the model which reduces to the Einstein-de 
Sitter model. The best radius of the boundary for fitting the m — z relation is found to be ~ 200 
Mpc. Moreover the latest data of a supernova with z ~ 1.7 (Riess et al. 2001) are naturally 
consistent with this last model with A = 0 (Tomita 2001b). For cosmological parameters in the 
inner region, we take the low-density value (Hq ~ 0.3), corresponding to the statistical analyses for 
local observations of galaxies and clusters. In the outer region, spatially flat models are chosen in 
accord with the recent results of CMB observations (Lange et al. 2001; Stomper et al. 2001; Pryke 
et al. 2001), and our models with small or vanishing A are found to be compatible with their data 
in the weak prior. 

Third, bulk flow on the scale more than 100 Mpc was found by by Hudson et al. (1999) and 
Willick (1999). It is impossible to account for its appearance as one of the linear perturbations in 
homogeneous cosmological models. In a model with a local void on scales ~ 200 Mpc, however, we 
can interpret the bulk flow as the nonlinear deviation of the local flow from the global expansion of 
the universe. This point is also one of the motivations for our inhomogeneous cosmological model. 

If a smooth power spectrum is assumed for the density perturbations, it is to be noticed here 
that the underdense nonlinear structure itself as a local void on scales ~ 200 Mpc is a unique 
object which evolved gravitationally from primordial adiabatic perturbations but whose amplitude 
deviates remarkably from that of average linear perturbations, so that it can be found at most only 
a few times within our horizon, and we cannot see anything similar in the neighborhood outside 
our local void. If similar structures were found more frequently, we would have a power spectrum 
with a sharp spike corresponding to their size. 

Owing to the gap of Hubble constants, an inner comoving observer has the additional velocity 
Vo [= 100(/i^ —/i^^)iio km/s] relative to the CMB (outer) frame, if the distance between the observer 
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and the center of the inner region is Rq Mpc. In a previous paper (Tomita 2000b) we derived the 
dipole moment D of CMB anisotropy caused by the additional velocity, and the Doppler velocity 
Vd corresponding to D, which depends on the model parameters in the inner and outer regions. 
If a certain condition for the parameters is satisfied, this velocity Vd gives the above-mentioned 
difference, so that we may have 

Vd = 166.6 Ud km s~^, (3) 

where Ud = (0.925, 0.124, 0.360) is also a directional unit vector. In §2 the expressions for 
the peculiar velocity in the linear approximation are derived in our model with a local void, and 
their difference from the homogeneous case is shown. In §3 the derivation of D and Vd is briefly 
accounted for and their numerical values are shown, and the above condition is discussed. Some 
mistakes found in a table in the previous paper are corrected there. In §4, it is shown that the bulk 
flows can be caused as the sum of the large-scale peculiar velocity and the additional velocity Vq- 
Section 5 is dedicated to concluding remarks. 


2. Peculiar velocity induced by inhomogeneities 

Our simple model with a local void consists of the inner homogeneous region (I^) and the outer 
homogeneous region (hii), whose boundary is spherical. The line elements are expressed as 

ds'^ = glMxTidx^r = -c\dRf + W{R)f{d{x^f + (4) 

where j (= I or II) represents the regions, f^{x^) = sinx^,X^ s-nd sinhy-^ for = 1,0,—1, 
respectively, and dQ'^ = dO^ + sin^ 9d(p^. The boundary shell is a time-like hypersurface S given as 
yi = Xi and = y^^^ 

The background models in and are specified by model parameters Hq, Dg and Ag 
for 3=1 and II, where Dg = [87rG/3(ilQ)^](/3o)'^ and Ag = and expressed using 

(= a^/(ao)^) and (= as 

dy^ /dR = {y^)~^/‘^Pj{y^), (5) 

where 

p,W) ^ K + W f + (1 - nj - (6) 

and Og is given by 

®o-^o ~ \/l — ^o~ ^0- 

Their Hubble constants and density parameters in the two regions are assumed to satisfy the 
relations Hq > and Dg < Dg^. 

In this section we consider local behaviors of matter motion in the Newtonian treatment. 
In terms of the common time coordinate t and Cartesian coordinates X, the equations of mass 



- 5 - 


continuity, motion and gravitation are 


dV 


dt 


|^ + Vjc(pV)=0, 

+ (F-Vjf)F + Vjf^> = 0, 


( 8 ) 


(9) 

^ = AttGp, (10) 

where p, V and ‘h are the mass density field, the velocity field and the gravitational potential, and 
X represents the proper distance. The background motion is expressed using uniform densities 
{p = p) and velocities {V oc X) in both regions. In order to describe local inhomogeneous motions, 
let us next use for X the comoving coordinates r in the inner region, where X = air and consider 
linear density and velocity perturbations (5i [= {p — pi)/pi\ and v in the inner region, where pi is 
the background (unperturbed) density in the inner region. Then the first two equations reduce to 


H-Vr • = 0, 

ot ai 

dv hi 1 „ , 

H- V H-Vr(/> = 0. 

at ai ai 


The gravitational equation is 


1 


Vr‘^(P = 4TTG{p- p), 


( 11 ) 

( 12 ) 

(13) 


which holds for p in both regions and the background density p is represented by p\.,pb and pn 
in the inner region, the shell, and the outer region, respectively. The solution for equation (13) is 
expressed as 


= -Gap 


i 


pi^i 


\r — r' 


r(fr' 


or 


= -Gai^ 




li \r - r' 


{^P)brb^ . , , 

— sm 0 do dp + 




'll k “ 


(14) 


(15) 


«I Jshell T - ’’b 

where Vb and {5p)b{^ {p — p)b) are the coordinate and density perturbation in the boundary shell 
respectively, and (h/o)ii = p — pii- 

Solving Eqs. (11) and (12), we obtain the velocity field in the inner region; 

G 


, , ai f 5i(r — r')d?r' 

v{r) = — + — 

4:71 J \r-r'\ «iJsheiiJo 


G 

+ - 

o-i Jii Jo 



piai^di 


^ Th^(v — 

ai{6p)bdt— -sin O'dO'dp' 


\r - rp 


i^p) 


II 


L i^p)i 


- 1 




\r — r 


/| 3 ' 


(16) 


where the first integral in the right-hand side of equation (16) is performed in both regions and we 
used the equation for density perturbation hi : 




( 17 ) 
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If we adopt the solution Di in the growing mode in equation (17), and use = dlnDi/dhiai, we 

obtain 


v{r) 


“I ttI fl 


6i{r — r')d^r' 


+ Res, 


r — r' 


(18) 


where 


Res = Ji 


I Jii \r-r'\'^ 


/shell \r - rb 


Ji 

J2 


G 


(«i) Jo 

G 

w/o 


ai{5p)idt, 


( (Mil _ \ 

V mi ) 


pia\ 5idt. 


(19) 


( 20 ) 


The first term (A) on the right-hand side of equation (18) represents the peculiar velocity induced 
in the model that is homogeneous in the entire region with and Hg. The second term {Res) gives 
the contribution of the local void to the peculiar velocity, which comes from the inhomogeneous 
matter distribution. In the homogeneous case without Res, the above expression accords with the 
usual one (e.g. Strauss & Willick (1995)). The peculiar velocity of the Local Group is given by 
the velocity at r = tlg as riG = 'I^(^lg)- 

The integrands in A and Res are proportional to \r — , and the main contributions in the 

integrations A and Res come from the region of distances ~ 40 Mpc and > 200 Mpc, respectively. 
Accordingly, Res may be smaller by the factor 5“^ than A, so that the rough estimate for the 
peculiar velocity is given by A and Res is negligible compared with A. 


Under the assumption that the small-scale density perturbations in the luminous matter is 
(3 X those in the dark matter and /3 is ~ 1, v{r) is calculated from equation (18) with Res = 0 
by investigating the inhomogeneity of spatial distributions in galactic redshift surveys. The values 
of i»(rLG) were derived by Rowan-Robinson et al. (2000) using the PSCz redshift survey including 
15,459 galaxies. When we average their values (in Table 1 of their paper) over the interval of 
|r| = (200 “ 300)h~^ Mpc, we obtain the expression in equation (2). 


3. Doppler velocity due to the CMB anisotropy for an off-center comoving observer 

In this section we consider the CMB anisotropy that is measured by an off-center observer 
in the inner region. For this purpose we must investigate the behavior of light rays in our model 
with a local void, by solving null-geodesic equations and extracting the multipole components. The 
detail is described in the previous paper (Tomita 2000b). Here we pay attention to light rays that 
are emitted at the recombination epoch in the outer region and reached the observer at the present 
epoch in the inner region. The redshift (^Jg^,) corresponding to the recombination epoch depends 
on the direction (j) between the light rays and 0(3. The value of zjgg is numerically calculated for 
Q < (f) < TT, and the dipole moment is derived. 
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When Zrec (= is given, the temperature T{(j)) of the cosmic background radiation is 

proportional to 1/(1 + ^^rec), and the dipole moment D and quadrupole moment Q are defined as 


077 /*27r 


D = 

= 2tt 


(l + Zrec) WlO sin (j)d(t>d(p\/{{1 + Zrec) ), 


0 ^0 


(1 + ^rec ) Wo sin (j)d(p\/{{1 + ^rec ) )? 


and 


077 r‘l77 


Q = 

= 27r 


(1 + ^rec ) sin(/d0(i(^|/((l + 

^rec) ) 


- 1 \ 


0 JO 


077 

/ (1 + Zrec)~^Y 20 sin (pdcj) /{(I + Zrec)”^), 
Jo 

where () means the average value taken over the whole sky, and 


^io(</>) = \l^cos4>, 


5 /3 



, - cos Q - 

47rV2 ^2 


Wo(0) = 

The Doppler velocity corresponding to D is given by 

Vd = c[{3/A7r)^^‘^D] 


( 21 ) 


( 22 ) 


(23) 


(24) 


This velocity Vd is a function of R (= UoXo)) i^o^ ^o) ^o^)- Here we fix as h} = 0.7, 

and fix Zb as Zb = 0.067, which corresponds to the case r;,(= a\xb) = 200(/i^)“^ Mpc. In the 6th and 
7th columns of Table 1, D and Vd for R = 10(/i^)“^ Mpc are displayed for various model parameters 
of the background models. First seven models are spatially flat in the outer region, and next three 
models are open. In the inner region all models are open. In the last column of Table 1 the values 
of R (= i7i7o) corresponding to Vd = 170 km s“^ are shown. It is found from this table that the 
distance R must be in the interval (10 - 20)(/i^)“^ Mpc. 

On the other hand, we assume that the Doppler velocity vector Vd gives the difference between 
ug and ucMB- From the comparison between uq in equation (2) and ucmb in equation (1), then 
the estimated Vd is given by equation (3). It is to be noticed that this velocity is in the direction 
of the Galactic center. 


The above derivation of D and Vd is quite the same as that in Tomita (2000b), but it was 
found that the expressions for D, Q and Vd in Table 1 of the previous paper have some mis¬ 
takes with respect to a factor 27r, and we should read D{x 10^), Q{ xlO^) and Vd (km s ^) as 
Zl(x 10^/27r), (5(x 10^/27r) and Vd (km s“^/27r) in Table 1 of Tomita (2000b). 


// Table 1 // 



4. Large-scale bulk flows 


At each point in the inner region we have a peculiar velocity v(r) given by equation (18) that 
is caused by the inhomogeneous matter distribution. Averaging the velocity over scales ~ Tbuik 
longer than 50/i“^ Mpc, we obtain the bulk flow velocity 'Ub(r*buik)- Branchini et al. (1999) derived 
the bulk flow velocity for IRAS galaxies in PSCz redshift survey including 11206 galaxies, measured 
out to 200h~^ Mpc. Its average value is 400 — 500 km s“^ for i?buik(= ao^’buik) = 60/i“^ Mpc and 
350 — 450 km s“^ for i?buik = 80h~^ Mpc. The amplitude and direction of bulk flows are sensitive 
to the measured scales of the matter distribution. If the scales are extended to more than 300/i“^ 
Mpc, the average values of the velocities may increase over the above ones. 

In our inhomogeneous models with a local void, we have not only the above bulk flow velocity 
Vb{r) corresponding to the peculiar velocity in the region V^, but also the velocity (ui) of the 
comoving observer (in V^) relative to the rest frame in which is the CMB frame. 

The latter velocity vi is equal to {Hq — H^)Rg in the direction CCi from the center to each 
galaxy G, where Rg is the distance between C and G. This velocity vi is expressed as the sum of 
the radial component (in the direction 0(i) and the constant component Vq (in the direction 
C(5), where |uo| = — H^)Ro and Rq is the distance between C and O, i.e., vi = Vy + Vq. The 

constancy of this component was shown in the previous paper (Tomita 2000b, 2002). The constant 
component Vq behaves as a peculiar velocity, while the radial component Vj. can be regarded as 
part of the Hubble motion from the viewpoint of the observer O. Accordingly, the total observed 
bulk flow velocity Vbuik is practically expressed using Ufe(rbuik) and Vq as 

"V^buik = ^’fe(r-buik) + Vo. (25) 

Here let us estimate this velocity, tentatively using one of the values of Branchini for i7f,(rbuik) on 
a scale of 60h~^ Mpc. Assuming the PSCz Ml method we adopt 

^’^(’’buik) = 450 rj-G km s“\ (26) 

whose direction is approximately taken to be equal to that of vq derived by Rowan-Robinson et 
al. (2000). We consider the case Hq = 0.3, = 1.0 and Aq = Aq^ = 0 and h^/hj) = 0.82, which 

is the best case for representing the SN data including the data of z = 1.7 (Tomita 2002). If we 
take Ro = 10//iq as in Table 1, we have 


Vo = 180 rid km s 

(27) 

Then we obtain from equation (25) 


^bulk — 378.3 TT-bulk km s , 

(28) 


where ribuik = (0.135, —0.906, 0.402) is a directional unit vector, so that the total bulk flow velocity 
has the amplitude 378 km s“^, which is 62% of ucmb, but its direction is comparatively close to 
that of ucMB in equation (1). 
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5. Concluding remarks 

In this paper we first derived the expression for peculiar velocities Uf,(r) in the inner region, 
which are caused gravitationally by linear matter density perturbations. Next we derived the 
Doppler velocity V(i{r) corresponding to the component of the CMB anisotropy caused by a local 
void. Using these results we obtained some constraints to our observer’s position O relative to the 
center C in the local void, under the assumption that Vd gives the difference between ucmb and 
ug (= 'I’(^’lg)) for the Local Group. It was found that the distance CO must be (10 - 20)/h} Mpc, 
depending on the parameters of the background models, and C(5 is in the direction of the Galactic 
center. 

Moreover, we estimated the amplitude and the direction of the bulk flow velocities on scales 
of ~ Mpc, as the sum of V}j{r) and the additional flow velocity velocity Vq in the direction 

of CCi which is specified so as to give the difference between the inner and outer Hubble constants. 

In our previous treatments (Tomita 2000b, 2002) we considered only the additional velocity 
to explain the origin of the bulk flow velocities obtained by Hudson et al. (1999) and Willick 
(1999). However, when we consider the peculiar velocity at the same time, the amplitude of the 
resultant bulk flow velocity is constrained by the observation so as to be much smaller than their 
flow velocities (~ 700 km/s), but the directions may be roughly consistent. 

In §5 we neglected Res for the estimation of vq, but because Res is of the order of 4% of uq, 
it may give to uq an error ~ 25 km s“^, so that it may bring an error ~ 14% to the difference of 
the velocities ~ 170 km s“^, or the distance GO (~ lQh~^ Mpc). 


This work was supported by Grant-in Aid for Scientific Research 12440063 from the Ministry 
of Education, Science, Sports and Gulture, Japan. 
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Fig. 1.— A simple model with a spherical boundary. A light ray is received by an observer O 
being near the center C, and Zb and z^ec are the redshifts at the boundary and at the recombination 
epoch, respectively, depending on angle (j). 
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Table 1; Dipole moment D and the velocity Vd in the case Zh = 0.067 and h} = 0.7. 



“0 


^0 


D (xlO^) 

Vd (km/sec) 

Rno{xh^) Mpi 

0.3 

1.0 

0.0 

0.0 

0.82 

11.76 

172.8 

bo 

0.4 

1.0 

0.0 

0.0 

0.82 

9.59 

141.0 

12.1^ 

0.5 

1.0 

0.0 

0.0 

0.82 

7.66 

112.6 

15.1^ 

0.6 

1.0 

0.0 

0.0 

0.82 

5.90 

86.7 

19.6^ 

0.3 

0.6 

0.205 

0.4 

0.80 

7.47 

109.9 

I 5 . 5 I 

0.3 

0.6 

0.303 

0.4 

0.87 

6.72 

98.9 

17.2^ 

0.4 

0.6 

0.269 

0.4 

0.82 

5.31 

78.0 

21.8^ 

0.3 

0.6 

0.0 

0.0 

0.80 

6.83 

100.4 

I 6.92 

0.3 

0.6 

0.0 

0.0 

0.87 

6.53 

95.9 

17.72 

0.4 

0.6 

0.0 

0.0 

0.82 

4.44 

65.2 

26.12 


^The outer space is flat. 
^The outer space is open. 



